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Abstract 

In three spatial dimensions, in the unitary limit of a non-relativistic quantum Bose 
or Fermi gas, the scattering length diverges. This occurs at a renormalization group fixed 
point, thus these systems present interesting examples of interacting scale-invariant models 
with dynamical exponent z = 2. We study this problem in two and three spatial dimensions 
using the S-matrix based approach to the thermodynamics we recently developed. It is 
well suited to the unitary limit where the S-matrix S = — 1, since it allows an expansion in 
the inverse coupling. We define a meaningful scale-invariant, unitary limit in two spatial 
dimensions, where again the scattering length diverges. In the two-dimensional case, the 
integral equation for the pseudo-energy becomes transcendentally algebraic, and we can 
easily compute the various universal scaling functions as a function of u/T, such as the 
energy per particle. The ratio of the shear viscosity to the entropy density n/s is above the 
conjectured lower bound of h/4irkB for all cases except attractive bosons. For attractive 2- 
component fermions, n/s > 6.07 h/Airks, whereas for attractive bosons n/s > QAh/AirkB- 
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I. INTRODUCTION 



In three spatial dimensions, in the unitary limit of a quantum Bose or Fermi gas 
with point-like interactions, the scattering length diverges. "Unitary" here refers to 
the limit on the cross section imposed by unitarity. These systems provide intriguing 
examples of interacting, scaling invariant theories with dynamical exponent z — 2, 
i.e. non-relativistic. The infinite scattering length occurs at a fixed point of the 
renormalization group in the zero temperature theory, thus the models are quantum 
critical. The only energy scales in the problem are the temperature and chemical 
potential, and thermodynamic properties are expected to reveal universal behavior. 
They can be realized experimentally by tuning the scattering length to ±00 using a 
Feshbach resonance. (See for instance and references therein.) They are also 

thought to occur at the surface of neutron stars. 

The systems have attracted much theoretical interest, and remain challenging 
problems due to the lack of small parameter for a perturbative expansion, such 
as n l / d a or n l / d r where a is the scattering length, and r the range of the 2-body 
potential. Early works were done by Leggett 

The universal scaling behavior was studied in {5, 6]. In 3 dimensions, this is the 
physics of the BCS/BEC crossover: since the fixed point occurs for an attractive 
coupling, the fermions may form a bosonic bound state which can subsequently 
undergo BEC. This cross-over was studied analytically in {^J by considering a model 



of 2-component fermions 



were performed in 
d < 4 



10, 
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coupled to the bosonic bound state. Monte-Carlo studies 



13J. The models can be studied in spatial dimension 2 < 



10| and an epsilon expansion carried out Il5j . There has also been some 



attempts to apply the AdS/CFT correspondence to these non-relativistic systems 
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19]. 

In the present work, we describe a new analytic approach to studying the unitary 
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limit based on our treatment of quantum gases in [20] , which appears to be well suited 
to the problem since it allows an expansion in the inverse coupling. Let us motivate 
this approach as follows. In one spatial dimension, the fixed point occurs for repulsive 



interactions. The model is integrable 21| and its thermodynamics determined exactly 



by the so-called thermodynamic Bethe ansatz (TBA) [22[ . In the TBA, the free energy 
is expressed in terms of a pseudo-energy which is a solution to an integral equation 
with a kernel related to the logarithm of the S-matrix. In the unitary limit the 
coupling goes to oo and the S-matrix S = — 1. The TBA is then identical to a free 
gas of fermionic particles. The formalism developed in {20 1 was modeled after the 
TBA: the free energy is expressed as a sum of diagrams where the vertices represent 
matrix elements of the logarithm of the (zero temperature) S-matrix. However since 
generally the N-body S-matrix does not factorize into 2-body scattering, the series 
cannot be summed exactly as in the TBA. Nevertheless, a consistent resummation of 
an infinite number of diagrams involving only 2-body scattering, the so-called foam 
diagrams, can serve as a useful approximation if the gas is not too dense. The result 
of summing these diagrams leads to an integral equation for a pseudo-energy, as in 
the TBA; in fact in 1 spatial dimension the TBA is recovered to lowest order in 
the kernel. Since the formalism is based on the S-matrix, it can be very useful for 
studying the unitary limit where S — — 1. 

In this paper we present the main formulas for the 3-dimensional case, however 
we mainly analyze the 2-dimensional case; analysis of the 3d case will be published 



separately 23j]. Phase transitions in two-dimensional Fermi gases were studied in e.g. 
251 ] . The fixed point separating the BEC and BCS regimes goes to zero coupling 
when d = 2, thus it is not obvious whether a unitary limit exists at strong coupling. 
As we will argue, there is a scale-invariant limit at infinite coupling g = ±00 where 
the S-matrix S = — 1. This is a meaningful unitary limit at very low (g = —00) or 
very high (g = +00) energy, although it does not correspond to a fixed point in the 



usual sense of a zero of the beta function. The scattering length indeed diverges in 
this limit. The possibility of this kind of unitary limit in two dimensions has not 
been considered before in the literature. 

In the next section we describe the unitary limit in 1,2 and 3 dimensions and 
its relation to the renormalization group. In section III, we define the interesting 
scaling functions for the free energy and single-particle energies by normalizing with 
respect to free theories. In section IV we describe the unitary limit of the formalism 



m [20( in both two and three dimensions, where the integral equation becomes scale 
invariant. The d = 2 case is especially simple since the kernel reduces to a constant 
and the integral equation becomes algebraic. Analysis of these equations in 2d is 
carried out for both infinitely repulsive or attractive fermions and bosons in sections 
V-VIII. The extension of our formalism to multiple species of possibly mixed bosonic 
and fermionic particles is considered in section IX. 

Kovtun et. al. conjectured that there is a universal lower bound to the ratio of 
the shear viscosity to the entropy density, 



where Jvb is Boltzmann's constant [24|. This was based on the AdS/CFT correspon- 
dence for relativistic theories in 3 spatial dimensions, and the bound is saturated 
for certain supersymmetric gau ge t heories. Counterexamples to the rj/s bound were 
suggested to be non-relativistic |24|, however no known fluid violates the bound. It 
thus interesting to study this ratio for non-relativistic theories, and in particular 
for 2d theories where no conjecture have been put forward. It has also been sug- 
gested that the gases in the unitary limit may represent the most perfect fluid, i.e. 
with the lowest value of r]/s. For the spacetimes considered thus far for a non- 
relativistic AdS/CFT correspondence in 2d, the result found is that rj/s is exactly 



1/Att 18|, |l9[ . We analyze rj/s for the attractive fermionic function of fi/T 



in section VI. If one disregards unphysical, potentially metastable regions, one finds 
rj/s > 6.07 h/4irkB- The other cases of repulsive bosons or fermions also satisfy the 
bound. However the attractive boson is below it: t]/s > 0Ah/4irkB- 



II. RENORMALIZATION GROUP AND THE UNITARY LIMIT 



The models considered in this paper are the simplest models of non-relativistic 
bosons and fermions with quartic interactions. The bosonic model is defined by the 
action for a complex scalar field 4>. 

S = J d d Kdt (#^0 - ^ - f (0W ) (2) 

(Throughout this paper, d refers to the number of spatial dimensions). For fermions, 
due to the fermionic statistics, one needs at least a 2-component field ^4.: 

S = J d d xdt ^ ~ - f$M^) ( 3 ) 

In both cases, positive g corresponds to repulsive interactions. 

The bosonic theory only has a U (1) symmetry. The fermionic theory on the other 



hand has the much larger SO (5) symmetry. This is evident from the work 26| which 



considered a relativistic version, since the same arguments apply to a non-relativistic 



kinetic term. This is also clear from the work 15] which considered an TV-component 
version with Sp(2N) symmetry, and noting that Sp(4) = SO (5). 

The free versions of the above theories have a scale invariance with dynamical 
exponent z = 2, i.e. are invariant under 

t ->■ hr 2 t, x ->■ A _1 x 

At a renormalization group fixed point, i.e. quantum critical point, the model is 
expected to have the same scale invariance. It is natural to define scaling dimensions 



dim[X] in units of inverse length or wave-vector k, i.e. dim[k] = 1, dim[x] = — 1, 
and dim[t] = —2. Requiring the action to have zero scaling dimension gives dim[0] = 
dimf^j — d/2 and dim[^] = 2 — d. The interaction is thus relevant for d < 2. 



A. 3d case 

The renormalization group behavior can be inferred from the coupling constant 
dependence of the S-matrix, or equivalently the kernels in section IV. For complete- 
ness we compute the beta function for arbitary d in Appendix A using conventional 
methods. Consider first a single boson in d — 3 dimensions. The kernel in eq. (145]) 
depends on the renormalized coupling gn given in eq. (I46p . where A is a high mo- 
mentum cutoff. The kernel G is related to the logarithm of the 2-body S-matrix as 
in eq. (HIT) , and from this we can deduce the S-matrix function: 

S(\k-k'\) = 167T/m9R -f - k '| (4) 
Vl w 167r/m£R + i|k-k'| v ; 

Galilean invariance, for equal mass particles, implies S is only a function of the 

relative momenta k, k' of the two in-coming particles. Unitarity of the S-matrix 

amounts to S*S = 1. Defining g = 'gA 2 ~ d , and requiring gn to be independent of A 

gives the beta-function: 

dg ~ m _ 2 

di = - 9 -4^ 9 (5) 
where I = — log A is the logarithm of a length scale. The above result agrees with 

the calculation in Appendix A. With the normalization of g in the fermionic case as 

given in ((3]), the resul ting beta function is the same as above. 

As discussed in (9), Il5| . the beta-function in d dimensions (see Appendix A) has 

the following implications. See Figure [TJ For d < 2, there is an infra-red stable fixed 

point at the positive value 



9* 



(2 - d)n d/2 T(d/2)2 d /m (6) 



corresponding to repulsive interactions. When d = 1, this fixed point is approached 
from g = ±00 when the cutoff A is infinite. In the limit of g — > ±00 the S-matrix, as 
given in eq. ((Hj) becomes S = — 1, which is the unitary limit. In the 1-dimensional 
thermodynamic Bethe ansatz, the thermodynamics is simply that of a free, one- 
component, fermionic gas, even though the original particle was a boson. 

For d > 2, there is an ultra-violet stable fixed point at the negative coupling g*, i.e. 
an attractive interaction. A bosonic gas with attractive interactions can be unstable. 
However for a fermionic gas, the attractive interaction can be balanced by the Fermi 
pressure, and the gas can be stable against collapse. 

We now describe the implications of the above fixed point on the scattering length. 
Consider a single boson. A straightforward calculation of the differential cross section 
in d spatial dimensions gives 

where ,M(|k — k'|) is the scattering amplitude. In the above formula k is the mo- 
mentum of one of the particles in the center of mass frame, i.e. 2k = |k — k'|. The 
Ai are an important ingredient of the kernels of the integral equations in section IV, 
and were computed to all orders in [20|] for instance. In d — 3, we equate a = no 2 . 
This gives 

a(k) = ™ - m (8) 
2vr ^1 + (mg R k/8n)* 



where g^ is the renormalized coupling If a(k) is measured at very small mo- 

mentum transfer |k — k'| ~ 0, this leads to the definition of the scattering length 

mgn = mg 
2tt 27r{l + mgA/47T 2 ) 1 ' 

where A is the ultra-violet cutoff. In order for a to diverge, g must be negative. In 
fact, at precisely the fixed point g* = g*/A, a — > ±00, depending on from which side 
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g* is approached. When g = g~, i.e. just less than g*, then a — > oo, whereas when 
g = g+, a — > — oo. The case g^ = — oo is on the BCS side of the crossover, whereas 
gn = +00 is on the BEC side. See Figured) 



d>2 



BEC BCS 



d=2 





.9* 



FIG. 1: Renormalization group flows; arrows indicate the flow to low energy. 

The S-matrix (j4j) has a pole at k = lQiri/mg^. Since physical bound states 
correspond to poles at \m(k) > 0, the bound state exists only on the BEC side of 
the critical point, with energy 

1287T 2 

-^bound-state — ^~ 5" (10) 

m?g l R 

In the BEC to BCS crossover from l/a = + to l/a = I/O - the binding energy goes 
to zero and the bound state disappears. Nevertheless, the crossover is expected to 
be smooth, if on the BEC side one includes the bound state in the thermodynamics. 



B. 2d case 



The d = 2 case is somewhat more subtle due to the marginality of the coupling 
g. The fixed point = and the RG flows are depicted in Figure [TJ The fixed 
point at <t* = is just the free field theory, thus there is no analog of the BEC/BCS 
cross-over. One can nevertheless formally define the unitary limit as S = —1, as in 
Id and 3d. In this subsection we explore this possibility and interpret it using the 
renormalization group. As we'll see, this limit occurs at g = ±oo, and the scattering 
length diverges. 

First begin with the beta-function in 2d, eq. (lllOj) : 

dg mg 2 



dlogA 4ir 

Let g = go at some arbitrary scale A . Integrating the beta function one finds: 



ii) 



^l-HfbgtA/Ao) (12) 
Note that g diverges at the scale: 

A, = A e 4w/m9 ° (13) 

This is the familiar Landau pole. Whereas in for example the relativistic 4 theory 
in d = 3 where the Landau pole is unphysical due to higher order corrections, here 
the beta function (11 II) is exact [15J, thus this divergence is physical. There are two 
cases to consider: 

Attractive case: If go is negative, A* < A , and thus g = — oo occurs in the infra-red. 
Repulsive case: If go is positive, A* > A , and g = +oo occurs in the ultra-violet. 
Both cases are consistent with the flows depicted in Figure [TJ 
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In fact, the scale A, = Ae 47r/m9 is an RG invariant: 

~n — r = ( 14 ) 
cHogA v ' 

and is the natural coupling constant in this problem that physical quantities are 

expressed in terms of. The S-matrix can be inferred from the kernel in eq. 

An/mg + log(2A/|k|) - in/2 



S(|k|) 

171 / Illy "|- lUg^ZiV/ T I'll / 4 

(15) 

lnaY9A /iWh - 7V/9 

(16) 



log(2A,/|k|) + z7r/2 

where |k| = |ki — k2 1 is the relative momentum of the two incoming particles. Near 
the fixed point g = 0~, A* = (A* = oo) and the S-matrix S — 1, consistent 
with the free theory. The same applies to the other side of the fixed point g = + 
where A* = oo. On the other hand, consider the low energy limit of the attractive 
case where g = — oo. Here A* = A, and this scale is effectively an infra-red cutoff. 
Thus at low energies | ki — k2 1 ~ 2 A* and S = —1. Similar arguments apply to the 
high energy limit of the repulsive case g = +oo. It is clear that, unlike in 3d, this 
definition of the unitary limit does not correspond to a renormalization group fixed 
point in the usual sense of a zero of the beta function and is somewhat engineered; 
nevertheless it defines a scale-invariant theory. 

Let us turn now to the scattering length. Using the scattering amplitude M. 



computed in 



20] , the cross section for relative momenta 2k is: 



Tn?g^ 1 

^ = ~Jk~ (l + IMl og{A/k)) 2 + 7l 2 ( 17 ) 

(18) 



Near the scale k = A*: 



k \og\Klk) + 167T 4 /mV 



2 /mgy 
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Equating a with a scattering length a, as appropriate to 2d, one sees that a diverges 
in the unitary limit g — > ±00. 

Because of the logarithmic dependence on | ki — k2 1 , the S-matrix does not have 
a pole. However the denominator is zero when |ki — k2 1 = 2k = 2iA*, which is a 
remanant of the bound state pole in 3d. The energy of this quasi-bound state is 
—2Al/m. Note that this quasi-bound state disappears as go — > 0~, since in this 
limit A* — > 0, and this is analagous to the situation in 3d. We thus expect that 
the "attractive" case of the unitary limit should be better behaved since this bound 
state disappears, in contrast to the repulsive case where A* — > 00 as go — > + , and 
this will be born out of our subsequent analysis. 

III. THERMODYNAMIC SCALING FUNCTIONS AT THE QUANTUM 
CRITICAL POINT 

At a quantum critical point, the only length scales of the quantum gas are the 
thermal wavelength = y^2n/mT and the length scale n 1 ^ set by the density 
n. Equivalently, one can express physical properties in terms of the only two energy 
scales, the temperature T and the chemical potential /z, since the density is a function 
of T, /i. 

In order to fix normalizations in a meaningful way, it is useful to consider the sim- 
plest theories with z = 2 scale invariance: free, non-relativistic bosons and fermions. 
We set k,B = fr = 1, except in a few select formulas. The free energy density is given 
by the well-known expression for each particle: 

^ = i/S log ( 1 - sc ~' J, " k ^ (20) 

where /3 = l/T, = k 2 /2m and s — 1,-1 corresponds to bosons, fermions respec- 
tively. The integrals over wave-vectors can be expressed in terms of polylogarithms 

11 



Li„(z), where z = is a fugacity, using J d d k = 2n d / 2 /T(d/2) J dkk d 1 and the 
integrals 



o 



00 zx u ~ l 

dx = T(v)IA v (z) 

e x — z 

00 7T v-i 

dx -r^- z = -r(^)Li,(-^) (21) 

valid for Re(z/) > 0. The result, as expected, is proportional to T/A^: 

fmT\ d/2 

r=-sT\—-\ U {d+2)/2 (sz) (22) 

There are two important limits to consider. For Bose gases near Bose-Einstein 
condensation, physically the interesting limit is fi/T — > 0. Since Li y (l) = C( z/ ) ) where 
C is Riemann's zeta function, this leads us to define the scaling functions Cd(fi/T): 

F=-^—cMT) (d=2) (23) 

C(5/2)m 3 / 2 T 5 / 2 
^ = (2^p ° 3(/i/ } ( = } 

where we have used £(2) = 7r 2 /6. With the above normalizations, Q = 1 for a single 
free boson when fi/T = 0. 

The above formulas are well-defined for fermions at zero chemical potential. Using 

-Li„(-1)= (l-^l)cH (24) 

one finds as fi/T — > 0: 

C2 = -, c 3 = 1 r= (free fermions) (25) 

2 2v 2 

It should be pointed out that the coefficients q are analogous to the Virasoro central 



charge for relativistic systems in d — 1, as discussed in [27]. 

The other interesting limit is T/fi — > 0, i.e. z — > 00 or z — > 0, depending on the 
sign of the chemical potential. Here the scaling forms are naturally based on the 
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zero temperature degenerate free fermion gas, where fi > is the Fermi energy. In 
fact, the function Li„(z) has a branch cut along the real axis from z = 1 to oo and 
the bosonic free energy is ill-defined at z = oo, in contrast with fermions. Using the 
analytic continuation of the asymptotic behavior 

- U^-z) « as oo (26) 

from positive integer v to half-integer values, we define the scaling functions bd as 
follows: 

T = b 2 (T/v) (d = 2) (27) 



4tt 

2V2m 3 / 2 ^ 2 
15vr 2 



The above normalizations are defined such that &<j(0) = 1 for a single free fermion. 
One can verify that the above d = 3 expression is the standard result for a zero 
temperature, single component, degenerate fermion gas where \x is the Fermi energy. 

Other thermodynamic quantities follow from the free energy. The pressure p = 
—T. The density is n = —dJ^/dfi and the entropy density is s = —dTj&T. Using 
the scaling form in eq. (1271) . one obtains for d = 3: 

2\/2 



ii 



157T 2 
2y/2 



(m/i)3/2 (l h - J h ') (28) 



157T 2 

where b' is the derivative with respect to its argument T/fi. (Henceforth, g' will 
always refer to the derivative of the function g with respect to it's argument fi/T or 
T/fi as defined above.) The analogous formulas in 2 dimensions are: 

n = ^ ( 2 b 2 - -b' 2 ) (29) 



Air \ fi 



s = -l—h 

47T 
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The energy density e follows from the relation e = Ts + \in + J 7 : 

6 = (30) 

It is interesting to note that the above result, which in terms of the pressure is 
simply e = pd/2, is a consequence of the mechanics of free gases. This shows that 
this relation continues to hold for interacting gases at a quantum critical point, as 
pointed out by Hoj^]. 

Also of interest is the energy per particle t\ = e/n: 

Consider first the limit T/fi — > 0, with [i positive. In this limit, if b is a smooth 
function of T/fi as T/fi — >■ + , then the b' terms in the density vanish and one simply 
obtains 

d 

lim ei = - — -fi (32) 
t/^o+ d + 2 v ' 

which is the same result as for a free gas, where for fermions /x is equal to the Fermi 
energy ep. The above formula is usually not appropriate to the T — > limit when 
/i is negative. For the interacting gas /i ^ €f, so this leads us to define the scaling 
functions £: 

um = d -^- 03) 

a ep 

As T/jU —7- + , the functions £<j should become universal constants, and for free 
fermions £d(0 + ) = 1. The Fermi energy ep can be defined based on its relation to 
density in the zero temperature free fermion gas. For bosons, one can formally use the 
same definition. For d — 2, ep — 2nn/m, whereas for d = 3, ep = (37r 2 n/ '\f2) 2 / z ' jm. 
This leads to the definitions: 

«■ - ^ < 34 > 

& = \ ("v) 5/2 Mr/ri 

ll 



Next consider the energy per particle in the limit fi/T — > 0. Here it is more 
appropriate to use the form in eq. fl23|) . which gives 

ti = ~T (35) 
z c d 

The expression for ei for free fermions in the limit fi/T — > in <i = 3 leads us now 
to define £: 

2 \2V2-2j C(3/2) WVW ; 1 ; 

With this normalization, for free fermions, as /x/T — >■ 0, £3 = 1. 

In two dimensions, e\ goes to zero for free bosons as fi/T — > since it is 
proportional to 1/^(1). However ei is finite for fermions in this limit. Using 
lim^ 2 (2 d/2 - 2)C(d/2) = 2 log 2, we define 

ei = ^6M (d = 2) (37) 

12 log 2 

With the above normalization, £2(0) = 1 for free fermions. 

Henceforth we drop the subscripts 2, 3 indicating the spatial dimension on the 
functions 02,03 etc., since in the sequel we will carry out the analysis of the d = 2 
case only. 



IV. THERMODYNAMICS FROM THE S-MATRIX IN THE UNITARY 
LIMIT FOR D=2,3 

In the formalism developed in [2^] , the contributions to the free energy can be ex- 
pressed as vacuum diagrams where the vertices are matrix elements of the logarithm 
of the S-matrix, and the lines are the filling fractions. There are vertices with 2N 
lines representing iV-body scattering for any N. The main result obtained in [20I . 
and reviewed below, is the formula for the free energy density in a self-consistent 
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resummation of the 2-body scattering "foam diagrams" . The 2-body scattering can 
be computed to all orders in the coupling in a standard zero-temperature calcula- 
tion. Clearly this result is still an approximation in that higher N-body scattering, 
and 2-body scattering contributions not of the foam diagram type, are neglected. 
Nevertheless, this approximation is expected to be valid when the gas is sufficiently 
dilute. 

The main ingredients of the formalism are as follows. Define the filling fractions 
as a function of a pseudo-energy e(k) 



which determine the density: 



d d k 1 

U ~ I (2^p e/M k ) - s (39) 



The consistent summation of 2-body scattering leads to an integral equation for the 
pseudo-energy e(k): 

e(k) = cu k - 11 - -log ( 1 + J (dk')G(k, k') e/fe(k ,)_ g j (40) 

The kernel G in this equation is related to the logarithm of the 2-body S-matrix S 
as follows: 

2it5(e-— (k 2 + k' 2 )^) VG(k,k') = ^(k.k'llog^^lk,^). (41) 
y 2m J 

where V is the spatial volume and (dk) = d d k/(27r) d . The kernel has the following 
structure: 

G = - j log(l + HM) (42) 

where Ai is the scattering amplitude and X represents the available phase space 
for two-body scattering. The argument of the log can be identified as the S-matrix 

16 



function. Finally, the free energy density then takes the simple form: 



^=-^ /W fslog(l + S /) 



2 Vl + s/o 

where /o is the filling fraction of the free theory: 

« k ) = JS^B _ , < 43 > 

Consider first the 1-dimensional bosonic case. The model is integrable, which 
implies that the N-body S-matrix factorizes into 2-body S-matrices, and the exact 
free energy is given by the thermodynamic Bethe ansatz 221 ] . The 2-body S-matrix 
is 

In the unitary limit g — > ±oo, S = —1, and the thermodynamic Bethe ansatz reduces 
to that of a free gas of fermionic particles. 

In 3 spatial dimensions the exact kernel for a single component boson is the 
following [201]: 



G(|k-k?|) 



1 - 


img R 
16tt 


k- 


k'| 


1 + 


img R 
16tt 


k- 


k'| 



m k — k' 



16tt 

arctan k — k 



(45) 



m|k — k'| V IQtt 

where gn is a renormalized coupling: 

1 1 mA , An . 

— = - + —. 46 

<7-R » 4vr 2 

with A an ultra-violet momentum cutoff. As in the 1-dimensional £lS discussed 
in section II, the unitary limit corresponds to gn ±oo where S — — 1. Thus in the 
unitary limit the kernel becomes 

G{k ' k ' ) = T m\k-k>\ {d = 3) (47) 
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where — , + corresponds to g being just below, above g*, where the scattering length 
a = +00 (BEC side) and a = —00 (BCS side) respectively. It should be kept in mind 
that the underlying interactions are attractive in both cases since the fixed point 
occurs at negative g. 



In two spatial dimensions the single boson kernel obtained in 20| is 



G(|k|) = log 1 



arctan 

m 



m M i+^(iog(t)+w2; 

mg/8 \ (48) 



arctan 



2tt 



m Vlog(2A*/|k|), 

where A* is defined in eq. (fT31) . and |k| = |ki — k 2 | is the relative momentum. In 
the unitary limit g — > ±00, the theory is at the scale A* and one should consider 
|ki — k 2 | ~ 2A + . The result is that G becomes a constant in this unitary limit: 

G(\k\) = T ^- (d=2) (49) 

lib 

In the attractive case, |k — k'| approaches 2 A* from above as g — > —00, and thus 
corresponds to the + sign above. The — sign then corresponds to the repulsive case 
where 2A* is approached from below. 

For two-component fermions, the phase space factors X in are doubled, and 
since G oc 1/X, the kernels have an extra 1/2 in the fermionic case: 

Cf crm i = -Gbose (50) 

The above unitary limit of the kernels leads to a scale-invariant integral equation 
for the pseudo-energy, which in turn leads to the scaling forms of the previous section. 
This will be described in detail for the d = 2 case in subsequent sections. Note 
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that the kernel has a well-defined expansion in 1/g; in the sequel we only consider 
the above leading terms. Solving the integral equation for intermediate values of g 
between and ±00 and inputting the solution into the expression for the free energy 
should track the RG flow between g = 0~ and —00, or between g = +00 and + , 
however we will not study this in the present work. 



V. ATTRACTIVE FERMIONS IN 2 DIMENSIONS 



The two-dimensional case is considerably simpler to analyze since the kernel G is 
a constant. In any dimension it is convenient to define: 

<f(k)=e(k) -£** + //; y(k) = e-^ (51) 

The expression ( l43|) for the free energy can be simplified to the following: 



-s log(l — se~ 



(52) 



2 e^-s 

Consider first a hypothetical single component fermion. In order to have a point- 
like local interaction one needs at least two components, and this will be treated in 
section VIII. As explained there, due to the SU(2) symmetry, the two-component 
fermion reduces to two identical copies of the following 1-component results. Since 
the kernel is independent if k, 5 is a constant. The integral equation then becomes 
the transcendental algebraic equation: 

y = 1 =f - log(l + zy) (53) 

y 

(We have used —Lii(—z) = log(l + z).) The — (+) sign corresponds to the repulsive 
(attractive) case. The free energy takes the form eq. ( 1231) where the scaling function, 
expressed in terms of the fugacity z, is 



~ (u 2 {-zy) + - y~ x ) log(l + zy)^j 



(54) 
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The scaling function b in eq. ([27)1 is 

3 log^ z 



TV 2 



b = T7T77T- C (55) 



and the density is 

TYlT 

n = — log(l + zy) (56) 
Zix 

The energy per particle scaling function £ is also expressed in terms of c: 

~ 2 log 2 
log(l + zy) 



Finally £ takes the form: 



log z 



2 



£ = r-rr t b (58) 

We will mainly plot these quantities against fi/T or its inverse; these functions are 

implicitly functions of the more physical quantity T/T F , where T F = ep/k B is the 

Fermi temperature. The relation with fi/T is 

T 1 

T F log(l + zy) 

There exists a solution to eq. (j55|) for any i.e. for all — oo < /i/T < +oo. The 
density is shown in Figure |2] and takes on all positive values and approaches oo as 
fi/T — > oo. As T//jL — > + , the density approaches the free field value n = m/x/27r. 
However at high temperatures there is a departure from the free field value: 

TflT 

lim n = 0.9546 (60) 

The filling fraction / is plotted in Figure |3j The scaling function c is shown in Figure 
HI One finds: 

lim c = 0.624816 (61) 

which is higher than the free field value c = 1/2. From these results, one obtains the 
equation of state: 

p= 1.077 nT, (p/T->0) (62) 
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One needs to make sure that all regions of fi/T are physical, and not for instance 
metastable. In particular, the entropy must increase with temperature, otherwise 
the specific heat is negative. In Figure [5] we plot the entropy density s, and one 
sees that ds/dT < when /i/T > 11.7. One must bear in mind that our formalism, 
though non-perturbative in some respects, is still an approximation, and such a 
feature could be an artifact that disappears once the corrections are incorporated. 
Since it is beyond the scope of this work to systematically explore these corrections, 
we will instead try and interpret our results as they are, for this and subsequent 
cases. The above feature of the entropy density could signify a phase transition at 
fi/T = 11.7, where the density and temperature are related by: 

TflT 

n = 13.1 T c = 0.076 T F (63) 

2n 

where Tp is the Fermi temperature ep/ks, where ep is defined in section III, and T c 
is the critical temperature for this hypothetic al p hase transition. The above value is 



comparable to the prediction T c « 0.1 Tp in [25] for quasi 2d (trapped) systems. In 
all other regions, and for all other cases considered below, ds/dT > 0. 

The scaling function b is shown in Figure El It has a well defined limit asT/fi — > 0, 
however the limiting value depends on whether \i is positive or negative: 

lim 6 = 0, lim b = 1 (64) 

Recall b = 1 is the free fermion value. The discontinuity at T j \i = is traced to 
the fact that c(oo) 7^ c(— 00), and has no physical significance. Since the density 
approaches m/x/27r, the equation of state is 

p = 27m 2 /m, {fi>/T — > 00) (65) 

The energy per particle scaling function £ is shown in Figured It has the limiting 
values: 

lim f= 0.907412, lim £ = 0.842766, lim £ = 00 (66) 
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Density for attractive fermions 

n/niT 




FIG. 2: Density of the attractive fermionic C3jSC clS cl function of \xjT . 

The above value for £ as /i/T — > — oo was determined numerically, however it turns 
out to equal 121og2/7r 2 . Equation ( 157)) then shows that the energy per particle 
ei = T, which simply means that in the limit \ijT — > — oo the gas is effectively 
classical. This feature will be encountered in other cases below. 

The single particle energy scaling function £ has the behaviour shown in Figure 
[SJ It has the limiting behavior 

lim £ = 1, lim f = 2.25593 (67) 

T/,u->-0+ T/n-^±oo 

The value £(0) = 1 is consistent with the arguments in [ji]]. It diverges as T//i-> CT 
due to the vanishing density at fi/T — >■ — oo. 
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P k 2 /2m 



FIG. 3: Filling fraction / for attractive fermions at two different values of fi/T. 
VI. RATIO OF SHEAR VISCOSITY TO ENTROPY DENSITY 



The ratio of the shear viscosity 77 to the entropy density s has units of h/ks in 
any dimension. A lower bound was conjectured for 3d relativistic systems 241] : 

h 



V 



> 

s Auks 



(68) 



based on the AdS/CFT correspondence. The bound is saturated in certain strongly 
coupled supersymmetric gauge theories. It is therefore interesting to study this ratio 
for non-relativistic, strongly interacting systems, in 3 and lower dimensions, since 
it is unknown whether there really is a lower bound. This ratio was studied for 3d 
unitary Fermi gases in 
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|-|32|. The attractive fermion case is the most interesting 
and well-behaved case in our formalism, as will be evident in the subsequent sections, 
so we study the issue in this case first. 

Consider first a gas with a single species of particle of mass m. For a non- 
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c for attractive fermions 

c 




FIG. 4: The scaling function c for attractive fermions as a function of fJ,/T. 
relativistic system in 2 spatial dimensions: 

V = ^rivm£ {ree (69) 

where v is the average speed and £f ree the mean free path. The mean free path is 
£{ ree = l/(\/2nc) where a is the total cross-section. The \/2 comes from the ratio 



of the mean speed to the mean relative speed[28(]. In the unitary limit, formally 
S — 1 + iXM. = —1, which implies the scattering amplitude A4. = 2i/I. (See 
section I V.) The cross-section in eq. ([7]) gives o = m 2 /kX 2 in 2 dimensions, where 



X = m/4 20j. Thus, in the unitary limit 
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° = Wi (70) 



where k is a single particle momentum. This gives 



= ( -mv 2 ) (71) 
1QV2 \2 
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Entropy for attractive fermions 
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FIG. 5: The entropy density s divided by m/j, for attractive fermions as a function of T/fj,. 

Since the pressure is determined by the average kinetic energy, and the ideal gas 
relation eq. ( 1301) still holds for a unitary gas, the average kinetic energy per particle 
is just t\ = —F/n, eq. (j3lj) . Finally we obtain: 

V rn FfdFY 1 (72) 



8 16V2n \3TJ 

where all the quantities in the above formula are the single component values of the 
last section. In terms of the scaling functions: 

rj 3 c n b 



s 4V27rC(2c-£c') 4V2b' {2b -lb 1 



(73) 



For two-component fermions the available phase space X is doubled. The cross- 
section is halved since spin up particles only scatter with spin down. Finally the 
entropy is doubled. This gives 

5) = 4 n (74) 

S ' fermi \ S / bose 
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FIG. 6: The scaling function b for attractive fermions as a function of T//I. 

The above expression is easily evaluated numerically using the expressions of the 
last section. The result is displayed for small values of fi/T in Figure [9j In this 
regime, rj/s is well above the conjectured bound, and comparable to the 3d values 
extracted from the experimental data 3]J] . We find 

lim 5 = 7.311-?- (75) 

fi/T->-0 S A-KKb 

This is well-below the values for common substances like water, however it is compa- 
rable to values for liquid helium, which is about 9 times the bound 24|. In the region 
shown: 

V - > 6.07-4- ( 76 ) 
s Attkb 

In Figure fTUl we plot rj/s as a function of T/Tp, and one observes a behavior quite 
similar both qualitatively and quantitatively to the 3d data summarized in 3l| . where 
the minimal value is about 6 times the bound. 
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Attractive^ f ermions 




FIG. 7: The energy per particle scaling function £ for attractive fermions as a function of 
»/T. 

Recall that in the last section, it was argued that the region fi/T > 11.7 is unphys- 
ical, or perhaps metastable, since the entropy increases with decreasing temperature 
there. Thus for the regions that are surely physical, the bound f!76|) holds. It is 
nevertheless interesting to study the viscosity in this unphysical region. In the zero 
temperature limit, i.e. Tj\i — > + , rj/s slowly decreases and seems to approach the 
bound. See Figure [TTJ However it eventually dips below it. See Figure O This 
behavior can be understood analytically as follows. For x = fi/T very large, the 
solution to the equation (15^]) is approximately: 

y(x) « y/x + 1/2 + (log x-l) fly/x (77) 

This leads to the asymptotic expansions: 

c w 3(x 2 + xlogx - x)/n 2 (78) 

TTlT 

n w — — (x + log(a;)/2) 

27T 
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Attractive fermions 




FIG. 8: The energy per particle scaling function £ for attractive fermions as a function of 

r//i. 



and 

7] 71 



S 2V2 

In terms of the density: 



(l/0og(Ai/r)-2) + r/2Ai) (79) 



lim ^ = ^= (80) 

m/t^oo s 2v^log(27rn/mT) 

It is important to note that although the energy per particle scaling function 
£ approaches the free field value at low temperatures, the above behavior is very 
different from the free field case. In the latter, the scaling function c = — 6Li 2 (— z)/tt 2 , 
which gives the diverging behavior at very low temperature: 

— ~ — -=- — , (free fermions) (81) 

s 4V27T T 

Finally, as \xjT — > —00, y ~ 1 + z, and c w Qz/n 2 and n ~ mTz/2ir. This leads 
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FIG. 9: The ratio rj/s as a function of fi/T. The horizontal line is l/4vr. 
to exponential growth: 

lim 0. = - jHL e -i>lT = _H^L I ( 8 2) 

m/t^-oo s 4 v / 2^1og(2vm/mT) 



VII. ATTRACTIVE BOSONS IN 2 DIMENSIONS 



Using the same definitions as for fermions in eq. (j51j) . the integral equation for 
bosons is 

y=l±-\og(l-zy) (83) 
V 

where the + (— ) sign refers to repulsive (attractive) interactions. The scaling function 
c is now 

c = A (hi 2 (zy) + ~(1 - y- 1 ) log(l - zyYj (84) 
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FIG. 10: The ratio r\j s as a function of T jTp for attractive fermions. 

The scaling function b has the same expression as in eq. fl55|) , with the above c. The 
density now is 

TTlT 

n = — — log(l - zy) (85) 
In 

The energy per particle scaling functions now take the form: 

~ 21og(2)c 
log(l - zy) 

and 

l0g ^ ^6 (87) 



log(l - zy) 

For this bosonic case, there is only a solution to eq. ( 183]) for z < z c ~ .34, or 
ji/T < —1.08. The density is shown in Figure [121 and note that it has a maximum. 
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FIG. 11: The ratio rj/s as a function of fi/T as /u/T gets very large.. The horizontal line 

is 1/47T. 

The filling fractions are shown for several fi/T up to logz c . From these plots, it 
appears that / is diverging at k = as z approaches z c , suggestive of condensation 
to a superfluid. Let us refer to the maximum density then as the critical density 

n c « 1.24^-^ <=^> T c «0.817> (88) 

The scaling function c is shown in Figure [151 The limiting value is 

c{z c ) « 0.35 (89) 

which is considerably less than for a free boson with c = 1. The function b decreases 
to zero at T/fi — > 0, and the limiting value is 

b(z c ) w 0.94 (90) 

The energy per particle scaling functions £ and £ are shown in Figures [Tj3, [T71 
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FIG. 12: The ratio rj/s as a function of f-t/T in the limit of very low temperatures (top 
curve). The horizontal line is 1/Att. The bottom curve is the approximation eq. (1791) . 



The limiting values are 



and 



lim f = 0.842766, £(z c ) « 0.39 (91) 

/.t/T— >— oo 



r lim _ £ = oo, f (* c ) ~ 2.3 (92) 



As for the fermionic case, the value £ = 0.842766 = 121og2/7r 2 implies the energy 
per particle e\ = T, which means the gas is in the classical limit. 

The ratio r]/s given by the expressions eq. ( I73p with the appropriate bosonic c 
and b. We find that r]/s has its minimum value at z c , where it is actually below the 
bound 

1 >0 A^— (93) 
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FIG. 13: Density of the attractive bosons as a function of \i/T . The origin of the axes is 
at (logz c ,0). 

VIII. REPULSIVE FERMIONS AND BOSONS IN 2 DIMENSIONS 
A. Fermions 

There are solutions to the eq. (|53|) for all fi/T < 0, and the density is positive 
in this range. A plot of the density is shown in Figure [THJ The density maximizes 
at fi/T PS —0.56, where n/mT ps 0.04. This seems physically reasonable given the 
strong repulsion in the unitary limit. In contrast, for small coupling g, the kernel 
G ~ — g, and there are solutions for positive chemical potential. In the bosonic case, 
our formalism leads to a critical density of n c = ^\og(27i/mg) for the Kosterlitz- 



Thouless transition 



20j. The filling fractions are shown in Figure [19j Note that 



they are considerably smaller than in the attractive expected. The scaling 
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f for attractive bosons 




FIG. 14: The filling fraction / for attractive bosons for fi/T = — 2, —1.5, —1.08 = logz c . 
function c is shown in Figure [20j It has the limiting value 

lim c = 0.303964, (94) 

which is considerably less than the free field value c = 1/2. The scaling function b 
has the limiting values lim T / M ^ - b = 0, and lim T / At ^„ 00 b = oo. 

Since the pressure is proportional to c, the same density occurs at two different 
pressures. One sees that for \ijT < —0.56, the density increases with pressure as it 
should. For fi/T > —0.56 the density decreases with increasing pressure, violating 
dp/dV < 0, and should thus be viewed as an unphysical region. We interpret this as 
a phase transition occuring at z c = e~°' 56 = 0.57, where the critical density 

n c = 0.25^fj ^ T c ^4.0 7> (95) 

At this critical point c(— 0.56) = 0.24. This perhaps corresponds to a transition to 
a crystaline phase. We cannot prove this, since we have not calculated the shear 
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FIG. 15: The scaling function c for attractive bosons as a function of fJ,/T. The origin of 
the axes is at (log .33, 0). 

modulus for instance, and it could simply be an artifact of our approximation, but 
let us take it as a hypothesis. Whereas a Wigner crystal phase occurs at low density 



in three dimensions, it occurs at high density in two dimensions [33|. The possibility 
of crystal phases for repulsive bosons was studied in {q]. For Coulomb repulsion with 
strength e 2 , the thermodynamic state of a classical Coulomb system is determined 
by T = \/rme 2 /kfiT, where n is the density Experimentally it was found that 
T PS 137[34j, which gives a critical density 

n c ^2.15x 10 9 T 2 ^— (96) 

On the other hand, since our model has point-like interactions, and is in the unitary 
limit, the conditions for a crystal phase are expected to be different. In particular our 
formula eq. ( 195]) has no e 2 dependence, and this leads to a linear in T dependence. 
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FIG. 16: The energy per particle scaling function £ for attractive bosons as a function of 
ixjT. 

For m equal to the electron mass 

n c « 4.5 x 10 9 T— ^- (97) 
cm 2 A 



It is interesting to note that for T of order 1, which is where the data in 34j was 
taken, the two densities ( 1961) and ( 1971) are comparable. 

The energy per particle function £ is shown in Figure [2TJ The limiting values are: 

lim £ = 0.842766, £(z e ) = 1.32 (98) 

fi/T->—oo 

As in previous cases, the above value of £ implies ei = T as fi/T — > — oo and the gas 
is thus in the classical limit. Beyond the critical point where the density increases, 
lim^/j^o <£ = oo. Note that £ starts to diverge around the proposed phase transition 
at fi/T = —0.56. Finally the other energy per particle function £ is less interesting, 
as it diverges at both endpoints of the density range where the density goes to zero. 
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FIG. 17: The energy per particle scaling function £ for attractive bosons as a function of 

r/ M . 

At the critical point it is still quite large: £(z c ) = 12.6. 

The ratio 77/5 reaches its minimum near z c , and is quite large compared to the 
attractive case: 

1 > 22.6 -V (99) 
s Ank B 

B. Repulsive Bosons 

The repulsive boson case is similar to the repulsive fermion, except that now 
there are only positive density solutions to the eq. fl83|) for /i/T < — log 2. Figure [221 
compares the density for bosons verses fermions. The maximum density for bosons, 
which could again possibly signify a critical point at z c ~ e -1,45 = 0.235, is half that 
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FIG. 18: Density of the repulsive fermionic function of \ijT . 

of the fermion case: 

n c = 0.125 <=^> T c ^8.0e F (100) 

At this point c(-1.3) « 0.11. 

The ratio rj/s is shown in Figure 1231 It has a minimum at fi/T = —1.72 near the 
critical point and in the physical region: 

^>8.85-V (101) 
s ~ Ank B 

Recall, the region ji/T > —1.45 is unphysical. 



IX. MULTIPLE SPECIES OF PARTICLES 

Let us suppose the gas consists of multiple species of particles labeled by the index 
"a", of mass m a , chemical potential /i a , and statistical parameter s a — ±1. Introduce 
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FIG. 19: Filling fraction / for repulsive fermions for fi/T = —3, —2, —0.55. The top curve 
corresponds to the critical density in eq. (|95|) . 



a pseudo-energy e a (k) for each type of particle and the filling fractions: 

fM = epJ) - s a (1 ° 2) 
As shown in Appendix B, in our approximation, the pseudo-energies satisfy the 
coupled integral equations: 

e (k) = w (k) - n a - - log (l + p G ab * J b b _ s ^ j (103) 

where w (k) = k 2 /2m a , 

Va = e~ PSa , S a = e a -u a + fi a (104) 

and G a b is the 2-body scattering kernel. By Galilean invariance, the kernel is a 
function of \v a — Vj,\ where v = k/m is a velocity. For any function g we defined: 

(G ab * g) (k) = J {db!)G ab % k')s(k') (105) 
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FIG. 20: The scaling function c for repulsive fermions as a function of n/T. 
The free energy density then takes the simple form: 



3a log (l 



-Pe a \ 1 Va 1 



') 



2(e^ - s fl 



(106) 



For the two-component fermion defined by the action ([3]), spin up particles scatter 
with spin down, and = G^. Thus, if the chemical potentials are equal, Ht = 
Hi = Hi then e-f = e±, and the thermodynamics is just two copies of the single 
component fermion described above. In particular, the density and c are doubled, 
but the energy per particle scaling functions are the same. 



X. CONCLUSIONS 



We have shown that the S-matrix based approach to quantum gases developed in 



20] leads to a new treatment of the scale-invariant unitary limit, where all of the 



thermodynamic scaling functions can be computed as a function of /x/T. Though our 
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FIG. 21: The energy per particle scaling function £ for repulsive fermions as a function of 
p/T. 

methods are of course not exact, they are sufficiently novel to provide new insights 
into these systems. It would be worthwhile to undertand the corrections to our results 
due to other types of diagrams involving N-body scattering for iV > 2, and also other 
more complicated diagrams involving 2-body scattering; this can be systematically 



explored using the full formalism in 



20]. 



In this paper we have mainly analyzed the 2-dimensional case, deferring the anal- 
ysis of the integral equations in the 3-dimensional case to a separate publication j^. 
For the 2-dimensional case, this required us to define a meaningful unitary limit 
where the S-matrix equals —1, and such a limit has not been considered before. 
We have calculated most of the interesting scaling functions for the free energy and 
energy per particle. 

The ratio of the shear viscosity to entropy density rj/s was also analyzed, and for 
fermions and repulsive bosons, it is above the conjectured lower bound of h/Anks- 

11 



Density for repulsive fermions verses bosons 

n/mT 




FIG. 22: Density of repulsive bosons (bottom curve) verses fermions as a function of n/T. 

For attractive bosons it drops below it by a factor of 0.4, however this could be an 
artifact of our approximation. For attractive fermions, the conjectured lower bound 
is reached at very large \ijT « 10 7 , however this was argued to occur in an unphysical 
or metastable region; in the physical regions, rj/s > 6.07 h/A-irks- 
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FIG. 23: The ratio rj/s as a function of fi/T for repulsive bosons. 
XII. APPENDIX A: THE BETA FUNCTION 

Consider a single free boson with action eq. §Z§. Setting the external energy and 
momentum to zero, the 4-point vertex at tree-level plus 1-loop is 

r« = - ig - tiEl I duddk ( • ^ ( * do?) 

y 2 J (27r) d + 1 \w-k 2 /2m + ie) \-uu -k 2 /2m + ie ) 1 ' 

The to integral can be done by closing the contour in the upper half plane picking 
up the pole at cu = — k 2 /2m + ie: 

mq 2 f d d k 1 
1 - loop = -i— / - — - — 108 
1 2 J (2n) d k 2 v ; 

Introducing an ultraviolet cut-off A: 

r « = -i(a m92Ad ~ 2 1 (109) 
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Defining the dimensionless coupling g as g = A 2 d cj, and requiring the I^ 4 ) be inde- 
pendent of A gives the beta function: 

where £ = — log A is the logarithm of a length scale; increasing I corresponds to 
flowing to lower energy. With the convention for the coupling g in eq. (J3]) for the 
fermionic case, the beta-function is the same as above but to half the symmetry 
factor of the 1-loop diagram. 



XIII. APPENDIX B: DERIVATION OF THE MULTI-COMPONENT CASE 



In this appendix, we extend the derivation in [20] to the general case of multiple 
species of particles, of possibly mixed Bose/Fermi statistics. Following the definitions 
in section VII, the 2-body foam diagram approximation is obtained by considering 
the free energy functional: 

r=-;/ W E(,.. 6 (i + a.)-^) (in) 

-iy(dk)(dk')J]7„(k)G oi ,(k,k')/ 1 (k') (112) 



where 



a.b 



The integral equation for the pseudo-energy e a follows from the variational principle 
6F/5f a = 0. Using 

|M = 5o , 6 (27r)^(k - kOr^-, (H4) 

OJb{ k ) 1 + s afo,a 

after some algebra one obtains the integral equation ( 11031) . (We have used 
G a b(k,k') = Gba(k',k).) Using the solution to the integral equation (II 031) in the 

44 



expression for F, one obtains the free energy density T = F in eq. (j!06p . 
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